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Abstract 

As is well known, the geometry of the interpolation site of a multivariate polynomial interpola- 
tion problem constitutes a dominant factor for the structures of the interpolation polynomials. 
Solving interpolation problems on interpolation sites with special geometries in theory may be 
a key step to the development of general multivariate interpolation theory. In this paper, we 
introduce a new type of 2-dimensional interpolation sites, tower interpolation sites, whose asso- 
ciated degree reducing Lagrange interpolation monomial and Newton bases w.r.t. fixed standard 
term orders such as lexicographical order, total degree lexicographical order, etc. can be figured 
out theoretically. Inputting these interpolation bases into Buchberger-Moller(BM) algorithm, we 
can also construct the reduced Grobner bases for related vanishing ideals. Experimental results 
show that in this way we can get the bases much faster than inputting tower sites directly into 
BM algorithm. 
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1. Introduction 

Let F be an arbitrary field and ¥ q a finite prime field of size q. IL d := F[xi, xn, . . . , Xd] 
stands for the d-variate polynomial ring over F. Given a set S = , . . . , ^ } C F d of 
[i distinct points. For prescribed values /,6F,! = l,...,/i, find all polynomials p G LT 
such that 

P{t®)=f u i = l,...,fi. (1) 
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The problem is called d-variate Lagrange interpolation, and the points are named inter- 
polation nodes, E the interpolation site, and p an interpolation polynomial. Generally, we 
only search interpolation polynomials in a subspace V C II . If there is a unique p £ P 
satisfying (1) for any fixed /i's, we call V an in terpolation spa ce that gives rise to a 
regular(ov poised) interpolation scheme (E,V)(see Lorentd (1993)). 

As a basic subject in Approximation Theory and Numerical Analysis, multivariate 
Lagrange interpolation has many applications in various fields of mathematics and com- 
puter science. Nevertheless, unlike univariate Lagrange interpolation, it is a relatively 
new topic with many basic problems unsolved. Since the set of all polynomials in Tl d 
vanishing at E constitutes the vanishing ideal 1(E) of E, multivariate L agrange interpo - 
lation becomes a particular exam ple of ideal interpolation introduced in Birkhofi ( 19791 ) 



and surveyed bv lde Boor J 20051). Recent year s , there has been considerable intere st in 
the topic, for example, cf. Sauer and Xu ( 1995 ); Bos et al. ( 2007 ); Shekhtman ( 20091 ) etc. 

In most cases, especially from a numerical point of view, we are not interested in all 
po lynomials sat isfying (1) but the one of "minimal degree", as in the univariate cases. 
In Sauer ( 20061 ). the notion of degree reducing interpolation spaces w.r.t. a fixed term 
order -< for interpolation problem (1) was introduced. Let ( E) C H d b e the subspace 
spanned by the Grobner escalier N^(I(S)) of 1(E) w.r.t. -<. Ide Boor! (|2007l ) pointed 



out that V^(E) is a canonical ^-degree reducing interpolation space as the unique one 
spanned by monomials. When d = 1, since 1(E) = ((x — £^)(x — £' 2 ') ■■ ■ (x — 
V^(E) = Span F {x fc : < k < fi - 1} is a classical one. But when d > 1, we can not 
obtain V^(E) theoretically for any E any more since the geometry of H, as a dominant 
factor for N^(I(S)), becomes very complex. In these cases, V^(E) can only be acquired 
by Buchberger-M611cr(BM for short) algorithm or variations thereof. 

As a mile stone algorithm for the c ompu tation of vanishing ideals, BM algorithm was 
presented in iMoller and Buchbergerl ( 1982 ) . For any interpolation site E and fixed term 
order -<, BM algorithm yields the reduced Grobner basis for 1(E) w.r.t. -< and the 
monomial basis(i.e. N^(X(5))) as well as a Newton basis for V^(E). However, as we all 
known, BM algorithm(and the variations) has a poor complexity that limits its appli- 
cations. One reason for this may be that the geometry of E is not concerned by these 
algorithms. Therefore, if we can solve the interpolation problems on some special dis- 
tributed nodes in theory, general interpolation theory will certainly benefit from it and 
the algorithms will also be impr o ved for some ca ses. 



In this decade, ISauerl (|2004 ); ICrainid (|2004) : IChen et al.l (|2006l) studied multivari 



ate Lagrange interpolation on a special type of interpolation sites, lower interpolation 
sites, and constructed the associated degree reducing interpolation monomial and New- 
ton bases in theory. Enlightened by them, this paper introduces a new type of bivariate 
interpolation sites, tower interpolation sites, whose geometry is more complex than lower 
sites nonetheless the related degree reducing interpolation bases can still be obtained 
theoretically. 

In this paper, we will study bivariate Lagrange interpolation on tower interpolation 
sites. After introducing tower sites and comparing them with lower sites in Section 3, we 
will present the degree reducing interpolation bases w.r.t. several standard term orders 
in Section 4 as main results. The last section, Section 5, will give our improved BM 
algorithm for tower sites and illustrate some experimental results. 

Next comes Section 2 that serves as a preparation for the paper. 
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2. Preliminaries 

For the reader's convenience, this section will introduce some notation and recall some 
basic facts. 

Let No stand for the monoid of nonnegative integers. A polynomial f G II 2 is of the 
form 



GN 2 :0^/ o - 6F} < 



The set of bivariate monomials in II 2 is denoted by T 2 . 

Fix a term order -< on II 2 that may be lexicographical order -<i cx , inverse lexico- 
graphical order -<i n i cx , total degree le xicographical order -<tdiex ; or to tal degree inverse 
lexicographical order -<tdinicx etc. (see lBecker and Weispfennind (1993)). For all nonzero 
/ G II 2 , we let LT(/) signify the leading term,LM(/) the leading monomial, and LC(/) 
the leading coefficient of /. Furthermore, for a non-empty subset Fell 2 , put 

LT(F) := |LT(/) : / G F}. 

Let A be a finite subset of N 2 ,. A is called a lower set if, for any G A, we always 
have 

{(*',/) GN 2 :0<i'<i,0<j'<i}c A. (2) 
Especially (0,0) G A. Let nij = maxr/jj)^ h, < j < z/ = max^ ^g^ fc. Obvi- 
ously .4 can be determined uniquely by the ordered (y + l)-tuple (mo, mi, . . . , m„) 
hence represented as L x (mo, mi, . . . , m„). In like manner, we can also represent A as 
h y (no, ni, . . . , n mo ) with n, = maxujA^ k, < i < mo- It is easy to see that v = uq. 



Definition 1 (|Crainid (|2004l )). We say that an interpolation site S in F is lower if there 
exists a lower set icN§ such that S can be written as 

Z = {(xi, yj ):(i,j)£A}, (3) 

where the x^s are distinct numbers, and similarly the y^s. We also say that S is .4-lower. 



In lde Boorl ( 2007 ). the -<— degree of a polynomial / G II 2 was defined to be the leading 
bidegrec w.r.t. -< 

S(f):=(i,j), xV' = LM(/), 

with i5(0) undefined. For any /, g G II 2 , if 6(f) -< (5(g) then we say that / is of Zower 
degree than <j and use the abbreviation 

f<9~ 5(f) ~< 6(g). 

In addition, / ^ g is interpreted as the degree of / is lower than or equal to that of g. 



Definition 2 ( Sauerl ( 20061 )). Fix term order -<. We call a subspacc V C II 2 a degree 
reducing interpolation space w.r.t. -< for the bivariate Lagrange interpolation (1) if 

(i) V is an interpolation space. 

(ii) V is -<— reducing, i.e., when L-p denotes the Lagrange projector with range V, then 
the interpolation polynomial 

L v q <q, Vqe II 2 . 
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If an interpolation space V for (1) is degree reducing w.r.t. -<, a basis {pi, . . . 
for V will be called a degree reducing interpolation basis w.r.t. -< for (1). Assume that 

Pi -< P2 -< < Pfi- If 

Pj(Z®) = Sy, l<l<j<fl, 

for some suitable reorder of S, then we call the basis a degree reducing interpolation 
Newton basis(DRINB) w.r.t. -< for (1). Let G^ be the reduced Grobner basis for the 
vanishing ideal w.r.t. -<. The set 

N„(X(3)) := {xV G T 2 : LT(g) \ x'y* , V 5 € G^} (4) 

is the Grobner escalier of 1(H) w.r.t. -<. Recall Section 1. As the monomial basis for 
■p^(H), N^(Z(S)) is named the degree reducing interpolation monomial 6asis (DRIMB) 
w.r.t. -< for (1). 

Finally, we will redescribe BM algorithm with the notation established above. 

Algorithm 1. (BM Algorithm) 

Input: An interpolation site S = {£W : i = 1, . . . , /i} c F d and a fixed term order -<. 

Output: The 3-tuple (G, N, Q), where G is the reduced Grobner basis for X(S) w.r.t. 
-<!, N is the Grobner escalier of X(H) (the DRIMB for (1) also) w.r.t. -<, and Q is a 
DRINB w.r.t. -< for (1). 

BM1. Start with lists G = [], N = [],Q = [],L = [1], and a matrix B = (6y) over F 
with pL columns and zero rows initially. 

BM2. If L = [], return (G,N,Q) and stop. Otherwise, choose the monomial t 
I'll / . and delete t from L. 

BM3. Compute the evaluation vector (£(£ < - 1 - ) ), . . . , t(£^)), and reduce it against the 
rows of B to obtain 

(«i, ...,«„) = (i(£ (1) ), . . . _ £ ai ( &ilj . . . , 6 V ), oi e F. 

i 

BM4.. If (^i, ... , v^) = (0, . . . , 0), then append the polynomial t — J^i a iQi to the list 
G, where qi is the ith clement of Q. Remove from L all the multiples of t. Continue with 
BM2. 

BM5. Otherwise (vi, . . . , u M ) ^ (0, . . . , 0), add (ui, . . . , v^) as a new row to B and 
t — aiqi as a new element to Q. Append the monomial t to N, and add to L those 
elements of {xii, . . . , Xdt} that are neither multiples of an element of L nor of LT(G). 
Continue with BM2. 



3. Tower interpolation sites 



Given an interpolation site E C F 2 . In ICrainid (|2004l ). two particular lower sets are 
constructed from H, denoted by S X {S), S v (z,), which reflect the geometry of H in certain 
sense. 

Concretely, we cover the nodes in S by lines 1$, If, . . . , I* parallel to the x-axis and 
assume that, without loss of generality, there are rrij + 1 nodes, say uf 7 -, 
on lj with mo > mi > • • • > m v 
Now. we define 



> hence the ordinates of ufj and ufj 



i ^ i', same. 



S X (S) := {(*, j) :0<i< mj , < j < u}, 



(5) 
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which clearly equals to h x (rriQ, mi, . . . , m v ). Similarly, we can also cover the nodes by 
lines ig, If, . . . , i^ parallel to the y-axis and denote the nodes on line l\ by u^ , u^, . . . , uf 
with no > n\ > • ■ • > n\ > hence the abscissae of and u^,,j ^ j' , same. Here, we 
put 

S y (E) := : < i < A, < j < n*} = L y (no,ni, . . . ,n A ). 

In addition, we can also define the sets of abscissae and ordinatcs 



Hj(S) :={x : (x, y) G if n H}, < j < i/, 
K(S):={y:(i,S)6lJ'nS} 1 0<*<A. 
Now, we will introduce our key notion: tower interpolation site. 



(0) 



Definition 3. We say that an interpolation site 3 in F 2 is x-tower if lower set S x (^) = 
La, (too, mi, . . . , m u ) C Nq with mo > mi > • • • > m v > such that 

E:={{ XiJ , yj ):(i,j)€S x (E)}, (7) 

where Xij £ Ho(E),(i,j) G 5 X (2), are distinct for fixed j. We also call S a ^(S)-^- 
tower interpolation site. In like fashion, if lower set S y (E) = L y (no,nx, . . . ,n\) C Ng, 
no > "-i > • • • > n\ > 0, such that 

S:={(x < ,v«):(i,j)G5 tf (S)}, (8) 

where G Vb(3), (i, j) G S V (E), arc distinct for fixed i, we will call S a y-tower inter- 
polation site. 

We have mentioned in Section 1, Definition 3 is enlightened by the notion of lower 
sites. Next, we will compare them in detail. Recall Definition 1. We have the following 
criterion. 

Theorem 4 ( Crainic ( 2004 )). An interpolation site 3 C F 2 is lower if and only if 
S x (E) = S y (E). 

Let 3 be an ,4-lower site. By Definition 1 and Theorem 4, we can deduce that A = 
S X (E) = S y (E). Hence, (3) can be redescribed as 

S-M:(<,j)6S.(S)} (9) 

or 

E:={{ Xi , yj ):(i,j)GS y (E)}. (10) 
Observing (7), (9) and (8), (10), we find that they are similar in form each to each 
nevertheless it seems that tower sites make the nodes on a line much freer than lower 
sites, which is supported by the following lemma as an alternative criterion for lower 
sites. 

Lemma 5. Resume the notation established above. An interpolation site 3 C F 2 is lower 
if and only if 

H (E)DH 1 (E)D---2H U (E) (11) 

or 

V (E)DV 1 (E)D---DV X (E). (12) 
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Proof. Assume that interpolation site S is lower. By (5), (6), and (9), we have Hj(E) = 
},j = 0,...,u. Since mo > mi > ■ ■ ■ > m v > 0, (11) follows. (12) can be 
proved in like manner. 

Conversely, we assume that (11) holds. Suppose S x (£) = L x (mo, mi, . . . , m v ) = 
Lj / (rig, n[, . . . , n'y) and S y (S) = L. y (rto, rii, . . . , n\). Since S X (S) is lower, we have A' = mo 
and n' Q = v. By (6), #-Hj(S) = m 3 + 1, j = 0, . . . , v. 

Next, we shall verify that A = mo = A', namely we can cover S by exactly mo + 1 
vertical lines l\ , i = 0, . . . , mo, which pass the mo + 1 nodes on line Iq respectively. Prove 
this by contradiction. Since A > mo apparently, equality can fail only when A > mo, i.e., 
there exists at least one node (x,y) £ H such that {x,y) ^ l\,i = 0, ...,mo- Because 
(x,y) £ 3, from (6), there must exist some -ffj«(5),0 < j* < v, such that x € /fj*(S). 
(11) implies that x £ H (E), which contradicts (5, y) £ V(,i = 0, . . . , m . 

Finally, we will prove 

n' h = n h , h = 0,...,m , (13) 
which implies ^(S) = S y (E>) immediately. We prove it by induction on h. When h = 0, 
for any node (x^,y^) £ l x v , x^ € ^(2) hence x^ € Hj{E),0 <j<v, due to (11). So 
we have found v+1 nodes on line x = x^ v ' that obviously equals to some l\* , < i* < mo. 
From the definition of v and S V (H), we have hq = n' = namely (13) is true for 
h = 0. Now, we assume (13) for < h < k < mo. When h = k + 1, by the induction 
hypothesis, we can find k + 1 lines, without loss of generality, that contain 

exactly rih + 1 nodes, < h < fc, of S respectively. Since iS x (S) = Lj / (n ,n' 1 , . . . ,nj„ ), 
we have (k + l,n' k+1 ) £ S X (E.) which implies that there must exist at least one node 
(x("'fc+i) ; yW+i)) e l x n , which is not on 1%, . . . , If. By (11), we can find n' k+1 + 1 nodes 

on line x = x^ rLk + 1 ' > that is a member of ■ ■ ■ , 'm } hence n^+i > n' k+1 . From the 

definition of 5* y (S), there exists one line in ■ ■ • ,lm a } that have exactly nk+i + 1 

nodes distributed among some lines of {Iq, . . . , l% k+1 } due to (11). Therefore, we have 
n 'k+i — n k+i that leads to rik+i = n' k+1 , namely (13) holds for h = k + 1. Consequently, 
Sa;(S) = L y (n , n' 1; . . . , «J no ) = L y (no, rii, . . . , Jt TOo ) = <Sy(5) which implies that 2 is lower 
due to Theorem 4. 

Swapping the roles of x and y, we can also prove the other statement similarly therefore 
complete the proof. □ 

Note that it is easy for Lemma 5 to extend to higher dimensions but difficult for 
Theorem 4. 

Now, we go back to tower sites. For an x-tower site S, although Hj(E) £, Ho(E),j = 
1, . . . , v, it may have nothing to do with any other i?i(S), i = 1, . . . , j — 1, j + 1, . . . , v. 
Comparing with (11), we find that, on one line, lower sites allow free amounts of nodes 
while tower sites allow free positions of nodes. The following proposition reveals the 
relation between them. 

Proposition 6. Resume the notation established above. 7/S satisfies 7?o(5) 2 ^i(^) 2 
■ • • 2 H„(S), then S is lower if and only if it is x-tower. Similarly, i/3 satisfies Vo(S) ^ 
Vi(3) 2 ' ' ' 2 ^a(S), then S is lower iff it is y -tower. 

Proof. Assume that 5 is lower with S x (£) = h x (mo,mi, . . . ,m v ). Since i?o(3) 2 
^i(3) 2 ' ' ' 2 H U {E), by (6), we have mo > mi > • • • > m„ > which implies that 



G 



S is a S x (S)-x-tower site according to (7). Conversely, if S is x-tower, since i?o(S) 2 
Hi(£) 2 ' ' ' 2 Hv(3), by Lemma 5, S is lower follows. 

In very like fashion, we can prove the other statement. □ 

Proposition 6 shows that the notions of lower sites and tower sites arc not mutually 
exclusive as there exist lower sites that are tower sites, but there are also tower sites are 
not lower and vice versa. The following examples show them with illustrations. 

Example 7. Observe (a), (b), (c) of Fig. 1. Si is an x-tower site that is not lower while 
S 2 is a lower site that is not x-towcr or ?/-tower. S 3 is a y-tower site that is lower too. 
(d) of Fig. 1 illustrates S y (S 3 ) that is equal to S^E^). 
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(c) 



(d) Sy(~ 3 ) 



Fig. 1. Illustrations for Example 7 
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4. Main Results 



In this section, we will pursue the degree reducing interpolation bases w.r.t. a fixed 
term order for bivariate Lagrange interpolation on a tower site E. Since an interpolation 
polynomial of minimal total degree is prefered in most situations, the -<tdiex and ^tdiniex 
cases will be discussed at the beginning. Note that we will continue with the notation 
established in the previous sections. 



Lemma 8 (jDong et al.l (|2005|)). Let E = {(xo, yo), ( x x,yo), (x m , yo)} CF 2 be a set 
of distinct nodes on line y = yo ■ Then 

1(E) = ((x - x )(x - xi) ■■■(x- x m ), y - y ). 

As a main result of the paper, the following theorem gives the DRIMB w.r.t. ^tdiex 
or -Ktdinlex f° r bivariate Lagrange interpolation on an x-tower or y-tower site respectively 
in theory. 

Theorem 9. Given an x-tower site S C F 2 . The DRIMB w.r.t. -<tdlex for (1) is 

N x = {x^ :(i,j)ES x (S)}. 
IfE is y-tower, then the DRIMB w.r.t. -< 

tdinlcx 

for (1) is 
Ny = {x i y* :(i,j)eS y (E)}. 

Proof. We will only prove the first statement. The second one can be verified in very 
like fashion. 

For fixed j, by Lemma 8, the ideal Xj := X({(x 0j ;, y 3 ), (sy, jjj), (x mj ,j,yj)}), j = 
0, . . . , equals 

((a; - x 0j )(x - xij) ■• • (x - x mji j),y- yj). 
Obviously, 2j's are pairwise comaximal hence the ideal 1(E) is equal to 

V V V 

j=0 j=0 j=0 

Let G be the reduced Grobner basis for 1(E) and Gk the reduced Grobner basis for 
[[3=0-^3 w - r -t- -^tdlexj < k < v. We will use induction on v to prove 

LT(G) = {x™ + \x m * +1 y,x™+ 1 1 ?,...,x m » +1 y",y»+ 1 }. (14) 
First of all, since 

2o = ((a; - x ,o)(x - x 1;0 ) ■ ■ ■ (x - x mo>0 ),y- y ), 

(14) is true for v = 0. 

Now, let v = 1. E is x-tower implies that mo > mi and H\(E) C Hq(E). Therefore, 

2b -2i =((x - x ,o) • • • (x - x mo>0 ),j/- y ) • ((x - x ,i) • • • (x - x mul ),y - yi) 
=((x - x , ) • ■ • (x - x mOi0 ), (x - x ,o) • • • (x - x mo , )(y - yi), 
(x - x ,i) • • • (x - x lnii i)(y - y ), (y - y )(y - yi)) 

=((x - x ,o) ■ ■ • (x - x m0i o), (x - xo,i) • • • (x - x mi ,i)(y - yo), (y - yo)(y - yx)}- 
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It is easy to check that LT(Gi) = {x mo+1 , x mi+1 y, y 2 } which means that (14) holds for 
v = 1. 

When v = 2, since mi > m 2 and i?2(S) C i3o(H), we get, after some easy computa- 
tions, 

2o -2i ■ ^2 =((x - xo,o) ■ •• (x - x TOo , ), (y - 2/0) (2/ - - 2/2), 

(x - X ,l) ■ ■ ■ (X - X mi> i)(x - £0,2) •••(£- 2Cma,2)(3/ - 2/o), 

(x - 2:0,1) • • • (x - x mi) i)(y - 2/0X2/ - 2/2), 
(x - 2:0,2) ■ • • (x - x m2 . 2 )(y - 2/0X2/ - 2/1)) 

/ (2) (2) (2) (2) (2), 

=:(.9o ,5i .92 .5s .94 >• 

Noticing that mi > 1112, we let q,f G IT 1 be the quotient and remainder respectively of 
the division of (x — xo.i) • ■ • (x — x mil i) by (x — xo. 2 ) • • • (x — x m2l 2), namely 

(x - x ,i) ■■■(x- x mit i) = q(x - 3:0,2) • ' • (x - x m2 , 2 ) + f. 

(2) (2) -(2) 

Denote the remainder of 173 w.r.t. g\ by #3 . One can check readily that 

-(2) (2) - (2) 

93 =93 -m ■ 



On the other hand, since 

(x - 2:0,2) ■ • ■ (x - x m2 ^)gf' - (x - x ,i) ■ • • (x - x mi .i)g^ 



9 2 



(2) _ L 



2/1 - 2/2 



we have g^ 9,3 94 — >+ 0. It follows that 

7-7-7- / _(2) „(2) -(2) (2)v 

Xo • Xi • X 2 — \ (7 '9l .93 .94 /• 

We claim that 

G 2 — G 2 .— |,g ,g 1 ,(,2/i-2/2j 93 .94 ). 

where 2/1 - 2/2 = LC(<^ 2) ), i.e., (2/1 - 2/2)~ 1 ffi is monic. 

In fact, if S(f,g) stands for the S-polynomial of arbitrary polynomials f,g £ II 2 , 

S(g^\g^) -— ^+ follows immediately because LM(g^) and LM(gp^) arc relatively 
prime. Notice that (x — 2:0,2) ■ • • (x — x m2 . 2 ) is a factor of (x — xo,o) ■ ■ • (x — x mo ,o), we get 

=2/ 2 . 9 ^ 2) - x m ™ 5 | 2) - (2/ - 2/ )(2/ - 2/i)9^ 2) + (y - Vo)(v - Vi)^ 
=9 2 9 2) - ^™9i 2) - (2/ - 2/o)(2/ - 2/1)9^ + ~ X - i 9l 2) 

(X - Xo,2) • ■ • (X - X m2 ,2j 
// , \ n (2) , ( (X - 2:0,0) • ■ • (X - X mo , ) mo _ m2 \ (2) 

=((2/o+ 2/1)2/ -yo2/i)9o + 7 \ } s-x - 94 • 

\{X - Xo,2) • • • (X - X TO2 ,2j / 

It is easy to check that LM(S(g { a 2) , gf ] )) = x m ° +1 y. Since 

LM(((2/ + 2/1)2/ - 2/o2/i)9^ 2) ) = z m ° +1 2/, 

LM ( ( (f-^0) -(f-^, ) _ x mo- m2 \ g (2)\ = x m 0y 2^ 
\\{X - Xo,2) ■ ■ ■ (X - X m2j2 ) / / 
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we have 

LM(S(gi 2 \ g| 2) )) = max (LM(((j/ + Vi)y - 2/o2/i).^ 2) ), 



-<tdlc 



(x - S , ) ■ ■ ■ (x - X mQ . ) _ ^ rlQ - m2 \ J2) 



LM f ^i22^ _ x™»-™= ^ 

V \(x - Zo,2j ' • ' [X - X m2t2 ) J 

=x mo+1 y : 

which implies that S(g^\ g^) — 0. Similarly, by 

S(g[ 2 \g {2) ) = (* m2+1 - (x - x 0>2 ) ... (i - . TO , 2 )) s f' - y 2 g| 2) , 

and LM(5(c?f \ 5 f ')) = x m * +1 y 2 , 

(LM((x m * +1 -(x- x Qft ) ■■■(x- * m2 , 2 )) 5l 2) ), LM(-y 2 g™)) 



max 

^tdlcx 



= max(x m2 y 3 ,a; m2+1 y 2 ) = x m2+1 y 2 = LM(S*(<?f \gf ] )) 

follows, namely S(gf\g^) 0. 
In like manner, we can also prove that 

S{g (2 \ fa - y^gf) 0, S(g[ 2 \ ( Vl y^g^) ^ 



0. 



Hence, there only remains S((y\ — y 2 ) 1 9^\g^) to be checked. Actually, it can 

be deduced that S((yi — y2)~ 1 g ( z\g^p) equals 

-t—gP + n-ii 

Vl - V2 



is x mi+1 y, therefore, from 

LM 





f (<7-x mi -™ 2 ) 5 f . 


-^ 2 U 2) ) 


, we know that its leadini 




= x m2 y 3 , 


- 2/2 / 






= x mi+1 y, 







S((yi-y2)- 1 gi 2 \g? ) )^+ 0. 

Ql 

The above arguments lead to the conclusion: S-polynomial S{g,g') — >+ for all 
g.g' G G 2 ,g ^ g' ■ Consequently, from Buchbcrger's S-pair criterion, G' 2 is a Grdbner 
basis for Iq ■ X\ ■ X 2 w.r.t. -<tdiex- 

Moreover, for any polynomial g £ G' 2 , it is evident that 

(1) LC(g) = 1, 

(2) No monomial of g lies in (LT(G 2 - {g})), 

which implies that G" 2 is the reduced Grdbner basis for Iq ■ Ty ■ X 2 w.r.t. ^tdicx- Hence 
(14) holds for v = 2. 
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Now, for v = k, assume (14). That is to say 

LT(G fe ) = {x m °+\ x m ^ +1 y, x m * +1 y 2 , x m » +1 y k , y k+1 }. 
Let G k = {g {k \... 

ifffc+i}- Without loss of generality, we may assume that 
LTfoW) =x m * +1 y\ i = 0,...,k, 

which imply that 

9o = {x- x 0:0 ) ■■■{x- x mo! o), 

g k % = (y - yo)(y - v\) ■ ■ • (y - Vk)- 

When v = k + 1, since S is x-tower, we obtain 
fc+i k 
J]l i =(IIW+ 1 

8=0 i=0 

= ( % ,---,9k+i) ■ ((x ~ xo,k+i) ■ ■ ■ (x - x mk+uk+1 ),y - y k+ i) 
l (fe+l) (k), \ 

g[ k '(x - £c ,A:+i) ■■■(x- x mk+lt k+i),g{ k '(y - yu+i), 

g { k-i( x - zo.fe+i) ■••(x- x mk+1 . k+1 ),g i k k } 1 (y - y k +i), 
9k ( x ~ x o,k+i) ■■■(x- x mk+1 . k+1 ),g[ k \y - y k +i), 
9k+i( x ~ Xo,k+i) ■•• (x - a:m fc+11 fc+x),5fc+i(j/ _ ffc+i)). 



where <7g + = <?o • Thus, <?q (y — y k +i) can be removed from the ideal basis above. By 
the induction hypothesis, we have 

si+t 1 ' := 9k+i(y ~ y fc +i) = (y ~ vo)(v -vi)---(y- 

We denote polynomial g^^x-xo^+x) ■ ■ ■ {x-x mk+1 ,k+i) by • For E x := {g[ k +i }, 

suppose g k k \y - y k +i) gt^ • Since m fe > m k+1 , we have 

LM(. 9 f +1) ) = x m " +1 y k - 
Furthermore, recalling the case v = 2, we can check easily that 

fffc fc) ( x ~ x 0,k+l) ■■■{X- Xm k+1 ,k+l) ^+ 0, 

where F\ := , 5j fc+1 ''}, which implies that the polynomial can be removed from 

the basis. 

Let E 2 := F\ and suppose g^l^y — J/fc+i) -^>+ ff^x • We can also deduce that 
LM( 5 fi + 1 1) ) = s^-i+iy*-! and 

ffi-ifa ~ »o,fc+x) ■•• (a; - Xm fe+1 ,fc+x) ^+ 0, 
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where F 2 := {g^\ g^ } . 

In this way we can construct two sequences E% , E2, . . . , -Bfc and -Pi , F2 , ■ ■ ■ , -Ffe such 
that 



(k), \ -Efc+i-i (fc+l) 

x - x 0>k+1 ) ■■■( 

where, for i = 1, . . . , k, 



li°\ x - Xo,k+i) ■ ■ ■ {x - x mk+1 ,k+i) Fk+1 '> +0, 



Fi - {9 k +2-V ■ ■ -:9k+l h 

*i-\9k+i ,9 k +i-il- 

Let <?i' c+1 ' ) = LC(gj fe+1 ) _1 g{ fe+1 . With quite similar methods that wc used in the case 
v = 2, we can prove finally that 

ll/O ' 9i 1 ■ ■ ■ i i/fc+i / 

is the reduced Grdbncr basis for Hi=o w.r.t. -<tdlex) Le., (14) holds when v = k + 1. 
By (4) and the definition of (S) , we have 

N^ tdlcx (J(S)) =N X = {aty ■■ G S,(S)}, 

which complete the proof. □ 

The next two theorems present degree reducing interpolation Newton bases w.r.t. 
^tdicx or -<tdinicx for bivariatc Lagrange interpolation on an x-tower or y-tower site 
respectively. 

Theorem 10. Let S C F 2 be an x-tower site of nodes u% ln = (x mn ,y n ), (m,n) € S X (B), 
which gives rise to polynomials 

3-1 t-l 
^ = <fij H(V- Vt) II^ - (hi) G S X (E), 

t=0 s=0 

where ipfj = l/II/Lilfe' — yt)Yi' s Ja( x ij ~ x sj) G F, and the empty products are taken as 
1. Then 

Qx = {<%:(?, j) eS w (E)} (15) 
is a DRINB w.r.t. ^tdiex for (1) satisfying 

<l>ij(Km) = S (i,j),(m,n), (hj) tinlex (m,7l). 

Proof. Fix (i,j) e S X (E). If (i,j) = (m,n), by y ^ j/i ^ • • ■ ^ % and x j ^ Zij ^ 
■ ■ ■ 7^ Xij, we have 

i-i i-i 

(=0 s=0 

Otherwise, if (i,j) >~inicx (wi, n), we have j > n, or j = n, i > m. When j > n, we have 

i— 1 

<%«wj = VijiVr, ~ Vo) ' • ' (j/n ~ ' • ' (j/n ~ Vj-l) J| (z mn ~ S s j) = 0, 

s=0 
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and when j = n,i > m, 

t=0 
n-1 

— J_ J_ ilJn Vtji^mn *^0n) ' ' ' (%vnn ^ran) ' ' ' i^mn %i— l,n) 

t=0 



0, 



which leads to 



<%«m) = °) (^i) >~iniox (m,n), 
namely Qj. is a Newton basis for Spanj-Q^.. By Theorem 9, it is easy to check that 
Span F (5x = SpanpATj; = V^ td]ox (3) . Therefore, Q x is a DRINB w.r.t. -<tdicx for (1). □ 

Similarly, we can prove the following theorem: 

Theorem 11. Let 3 C F 2 be a y-tower site of nodes u v mn = (x m , y mn ), (m, n) G S y (E). 
We define the polynomials 

i-i 3-1 

s=0 t=0 

where ip^ — 1/ YY s Jo^ Xi ~ x s)Y^t=o{Vij ~ Da) S F. The empty products are taken as 1. 
Then, 

Q y = {4>l:{i,j)eS v {E)} (16) 

is a DRINB w.r.t. -<tdinicx for (1) satisfying 

<t>ij(Umn) = %,i),(m,n), (hj) ^icx (m,n). 

Now, we will turn to -<i ox and -<i n icx cases. For any interpolation site 3 C F 2 , the 
following theorem presents the degree reducing interpolation monomial and Newton bases 
w.r.t. ^icx or -<i n iex for (1) theoretically. 

Theorem 12 (jWang et al.l (|2009l )). Let 

s = {Km = Kn.J/D : ( m , n ) & S X (E)} = {u v mn = (x v mn , y^J : (m, n) e S y (E)} 

be an interpolation site in F 2 . Then 

(i) the set N x := {x % y' : G S X (E)} is the DRIMB as well as Q x := {(% : G 
S^S)} is a DRINB w.r.t. ^i cx for (1), where 

3-1 i-l 

ffij = II (f ~ Vat) II ( x " <J )' (*' J) e S * ( S ) ' 

t=0 s=0 

with ipfj = l/Y\t=o(Voj — Vot) YVs=o( x ij ~ e ^ ^ e em pty products taken as 1; 

(ii) tfte set N y := {xV : G S y (E)} is the DRIMB as well as Q y := {0V. : g 
5y(5)} is a DRINB w.r.t. -<i n icx /or (1), where 



i-l 3-1 



s=0 t=0 
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with cpfj = 1 / YVs=o( x io ~ x so) Ylt=o(yij ~ Vit) an< ^ the empty products are taken as 
1. 

Since Theorem 12 holds for any interpolation site in F 2 , it obviously holds for tower 
sites. It is easy to find that the N x in the theorem is same to the one in Theorem 9, 
namely the DRIMB w.r.t. -<x ex is same to the one w.r.t. -<tdiex for interpolation on an 
x-tower site. In addition, the Q x here is same to (15) when 2 is x-tower. For N y and Q yi 
similar statements are also true for a y-tower site. 

In a word, for a preassigned x-tower site, we can obtain the DRIMB w.r.t. -<tdicx or 
^iex by Theorem 9, and a DRINB w.r.t. -<tdiex ° r ~<icx by Theorem 10. The y-tower cases 
can also be solved by Theorem 9 and 11. 



5. Algorithm and Experimental results 

Let 5 be a tower site in F 2 . Suppose (G, N, Q) be the 3-tuple output of BM algorithm. 
Theorem 9-11 prese nt us N and Q theor etically hence we can obtain them with little 
effort. According to Marinari et al. ( 19931 ). LT(G) is contained in the border set of N. 
Therefore, inputting N, Q into BM algorithm can make us obtain G much faster than 
compute G directly with BM algorithm. 

Algorithm 2. (TBM algorithm) 

Input: A x-towcr(y-tower) interpolation site S C F 2 of \i nodes and fixed term order 

~<tdlex(~<tdinlex) Or -<lex(^Iinlex)- 

Output: The 3-tuple (G, N, Q), where G is the reduced Grobner basis for X(S), N is 
the Grobner escalier N(J(E)), and Q is a DRINB for (1). 

TBM1. Construct lower set S x (E)(S y (S)) following the process introduced in Section 



TBM2. Compute the sets N and Q according to Theorem 9-11. 

TBM3. Compute the border set L := {x ■ t : t e N} [j{y ■ t : t e N} \ N . 

TBM4. Construct fix /i matrix B whose (h, k) entry is (uf ){<f>% (uT)) where </>^(0^), 
u%(u v k ) are hth and fcth elements of Q = {4>ij(<t>ij) '■ {hi) € S x (E)(S y (3))} and 3 = 
{u^Ju^J : (m,n) g S x (E)(S y (E))} w.r.t. the increasing -<miex(-<lex) on and 
(m, n) respectively. 

TBM5. Goto BM2 with N, Q, L, B for the reduced Grobner basis G. 

In the following, we will show the timings for the computations of BM-problems on 
tower sites in finite prime fields w.r.t. term order -<inicx or -<tdiex- BM and TBM algo- 
rithms were implemented on Maple 12 installed on a laptop with 2 Gb RAM and 1.8 
GHz CPU. 

For field F 23 and -< t diex, 



Algorithms^ 


#3 


100 


300 


500 


900 


1300 




TBM 




0.125 s 


1.297 s 


4.953 s 


20.609 s 


56.906 s 


BM 




0.531 s 


9.437 s 


53.625 s 


224.812 s 


861.531 s 



For field F43 and -< 



MnlCX; 



14 



Algorithms 


#2 


600 


900 


1200 


1500 


1800 


TBM 




7.781 s 


21.094 s 


45.282 s 


87.734 s 


134.265 s 


BM 




121.094 s 


362.000 s 


1002.578 s 


2077.907 


2729.750 s 
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